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Abstract 

We present a direct derivation of the thermodynamic integral equations of the 0(3) 



nonlinear cr-model in two dimensions. 
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The two-dimensional 0(3) nonlinear cx-model is a well known example of an exactly 
integrable field theory |3| which is asymptotically free at short distances and a mass gap 
appears in the spectrum through the dimensional transmutation mechanism Q. Its action 
A is defined by 

A—-[ d 2 x(d u ff) 2 (1) 
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where fj = (n x , n y , n z ) is a three component unit vector on the sphere S 2 ; ff = 1. 

It is believed that the 0(3) nonlinear u-model has no bound states and its fundamental 
particles belong to the 0(3) multiplet. Using these hypotheses and the factorized bootstrap 
A. Zamolodchikov and Al. Zamolodchikov have found the respective two-body S-matrices 

smpy, 



The solution for the 0(3) nonlinear cr-model directly from its bosonic representation (1) 
is still unknown. However, Wiegmann ||, using a fermionic formulation proposed early by 
Polyakov and Wiegmann || , was able to find the associated Bethe ansatz equations from a 
subtle limit of an exactly solved four fermion model. These equations are given by, 

e*"^ = ft ^~^ + m \[ ^~ X ' + m i < = 1,2,...,JV (3) 
j= \;^ A - Pj - in £=1 A - Aj- - itt 

n ^ r:: u j ; f ■ , = i,2,..,m (4) 

i= l,7^ A* - Aj- - Z7T /?i - Ai - 27T 

where L is the size of the system, m is the mass of the particles A^i = 1,2,3 ; and the 
corresponding energy is E = J2{Li mch{0i). 

Recently there has been a renewed interest in studying several properties of model (1). 
M. Liischer and U. Wolff have measured the S-matrices (2) , exploring the two-particle 
phase-shift. P. Hasenfratz et all and P. Hasenfratz and F. Wiedermayer |J have calculated 
the exact value of the mass gap m in terms of the renormalization group scale. V. Fateev and 
Al. Zamolodchikov have conjectured the thermodynamic Bethe ansatz (TBA) equations 
of the nonlinear 0(3) a-model as the M — ► oo limit from those of the Z(M) parafermionic 
field theory perturbed by a certain relevant field that breaks the Z(M) symmetry. 
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In this letter we present a direct derivation of the TBA equations of the 0(3) nonlinear 
cr-model from Eq.(2). This is a reliable approach and has the advantage of allowing one to 
generalize our results in the general case of 0(N) symmetry. We also discuss the presence 
of a magnetic field coupled with its conjugate conserved Noether charge. 

In order to deduce such TBA equations, one has first to notice that the eigenvalue 
problem for the S-matrices (2) defined by finding a matrix R(A) satisfying the relation, 

R(X - u)T(X) (g) TO) = TO) (g) T(X)R(X - u) (5) 

where the monodromy matrix T(A) is 

T(A) = S{p x - X)S(p 2 - X)....S(p n - A) (6) 

is related with the similar problem solved in the case of the Heisenberg chains [|10||. Indeed 



the rapidities A-,-, satisfying Eq.(4), can be interpreted as the massless excitations of the spin- 
1 Heisenberg chain [110] . This observation allows us to conclude that in the thermodymanic 



limit, L — * oo, the rapidities Xj form strings of the following type [11]: 



ITT 

\ j = XJ 1 + —(m + 1 - 2a); a = 1,2, ...,m (7) 

where the AJ 1 are real; the center of the m-string Xj. 

Substituting Eq.(7) in Eqs.(3,4), taking the logarithm and the thermodynamic limit , 
we rewrite Eqs.(3,4) as, 

/OO °° POO 

* X (J3 - f3')a((3')df3' - £ / - X)p(X)dX (8) 

-oo m=l J -°° 

00 POO POO 

27rJ m (A) + ]T / Q mk (X - X') Pk (X')dX' = / * w (/3 - X)a((3)d(3 (9) 

k—l J —oo J —oo 

where ™ = a{(3) + a((3), ^ = Pm {X) + p(X); a{(3) {a{{3)) and p(A) (p(A)) are the 
densities of particles (holes) related with the rapidities (3j and Xj, respectively. The functions 
^m(x) and Qmk(x) are defined by, 

*m0) = 9m+i(x) + g m -i{x) (10) 



©mfcO) = J2(gi(x) + g i+ 2{x))\ i= \m-k\,\m-k\ + 2,...,m + k-2 (11) 

i 

where g a (x) = 2artg(^) and g Q {x) = 0. 

The temperature is encoded by using the method pioneered by Yang and Yang in the 
case of the delta function interaction ||12|| . This technique consists in the minimization of 
the free energy F = E — TS in the equilibrium state. In our model, the energy E and the 
entropy S are given in terms of the densities a((3),a(/3), p(X), p(X) by, 

/oo 
mch(x)po(x)dx (12) 
-oo 



S = [(Pm(x) + p~ m (x))ln(p m (x) + p m {x)) - p m (x)ln(p m (x)) - p~ m (x)ln(p~ m (x))) 

(13) 

where we made use of a simplified notation, a(x) (&(x)) = po(x) (po(x)). 
Using 5F = 0, we find the following equations 

m 00 
X (f3) = -ch((3) - - (3') * ln(l + e -*P)) - £ K(P ~ A) * ln(l + e~^)(U) 

1 m=l 

oo 

ln{l + e em(A) ) = V m (\ - p) * ln(l + e"*^) + ]T A km (\ - A') * ln(l + e -^) (15) 

k=l 

where we have redefined the densities ? = e _x ; ^ = e~ 6m . The symbol g(x — y) * f(y) 
denote the convolution ^ /f^ #(x - y)f(y)dy\ W m (x) = and A mfc (x) = £<d m k(x) + 

2n5 mik 5(x). 

In order to simplify the Eqs.(14,15) it is important to use the following relations, 

A mk( x ) = 27r5 m>fc 5(x) - (f)(x)(5 m+1>k + S m -i )k ) (16) 

*' m (x) = c/)(x - y) * A m2 (y); = ^-y (17) 

First applying A km in Eq.(15) and eliminating the sum in Eq.(14) by using Eq.(15), we 
finally obtain the following result, 

X {fi) + 4>{P - A) * ln(l + e"«W) = jch(P) (18) 



e m (A) + (f>(\ - A') * [ln(l + e^ 1 ^) + ln(l + e'^- 1 ^)] + 5 mj2 (f>(\ - (3) * ln{\ + e~ xW ) = 

(19) 

in accordance with the results conjectured in ref. ||. 

It is also interesting to study the model (1) in presence of a magnetic field h coupled with 
one of the components of the conserved charge Q = J fjXd t fjd 2 x. Now, the energy changes 
to E = J2f=i mch(Pj) + h(N — M) and the minimization procedure follows close to the one 
for the Heisenberg chains [fTC] , |10| . In this case, we have an extra boundary condition to the 
variables e m , such that linim^oo ^ = h. Strictly at T=0, the rapidities Xj are frozen, and 
one can reproduce the known results obtained by Wiegmann 0. 

It will be interesting to generalize our results to the general O(N) cr-model. In doing 
that, we might generate some new TBA equations for perturbed conformal field theories, 
by projecting out the strings at some finite value m. Such approach has successfully been 
applied in the case of the critical RSOS model |0|. Based on ref. @ and our present result 
this approach can also be used in the general O(N) a-models. 
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